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In the atmospheric plasma of a strongly magnetized neutron star, vacuum polarization can induce 
a MSW-like resonance across which a X-ray photon may (depending on its energy) convert from 
one mode into the other, with significant changes in opacities and polarizations. We show that 
this vacuum resonance effect gives rise to an unique energy-dependent polarization signature in the 
surface emission from neutron stars: for "normal" field strengths (10 12 < B < 7 x 10 13 G), the 
plane of linear polarization at photon energy E < 1 keV is perpendicular to that at E > 4 keV, 
while for "superstrong" field strengths (B > 7x 10 13 G), the polarization planes at different energies 
coincide. The detection of polarized X-rays from neutron stars can provide a direct probe of strong- 
field quantum electrodynamics and constrain the neutron star magnetic field and geometry. 



Surface emission from isolated neutron stars (radio pul- 
sars and radio quiet neutron stars, including magnetars 
0) provides an useful probe of the neutron star (NS) in- 
terior physics, surface magnetic fields, and composition. 
The advent of X-ray telescopes in recent years has made 
the detection and detailed study of NS surface emission 
a reality 2] . On the other hand, X-ray emission from ac- 
creting X-ray pulsars has long yielded information on the 
dynamics of accretion and radiative mechanisms in the 
strong gravity, strong magnetic field regime 0]. 

In the magnetized plasma that characterizes NS atmo- 
spheres, X-ray photons propagate in two normal modes: 
the ordinary mode (O-mode) is mostly polarized par- 
allel to the k-B plane, while the extraordinary mode 
(X-mode) is mostly polarized perpendicular to the k- 
B plane, where k is the photon wave vector and B is 
the external magnetic field [3j . This description of nor- 
mal modes applies under typical conditions, when the 
photon energy E is much less than the electron cy- 
clotron energy Es e = TieB / (m e c) = 11.6 B\i keV [where 
B12 = B/(10 12 G)], E is not too close to the ion cyclotron 
energy Eb% = 6.3 B\2(Z/A) eV (where Z and A are the 
charge number and mass number of the ion) , the plasma 
density is not too close to the vacuum resonance (see be- 
low) and OkB (the angle between k and B) is not close to 
zero. Under these conditions, the X-mode opacity (due 
to scattering and absorption) is greatly suppressed com- 
pared to the O-mode opacity, kx ~ (E / Ese) 2 k<j 0,13- 
As a result, the X-mode photons escape from deeper, hot- 
ter layers of the NS atmosphere than the O-mode pho- 
tons, and the emergent radiation is linearly polarized to a 
high degree (as high as 100%) [HIHOl- Measurements of 
X-ray polarization, particularly when phase-resolved and 
measured in different energy bands, could provide unique 
constraints on the magnetic field strength and geometry 
and the compactness of the NS |E S U- 

It has long been predicted from quantum electrody- 
namics (QED) that in a strong magnetic field the vac- 
uum becomes birefringent 0, [lJJ, 0, E3- Acting 
by itself, the birefringence from vacuum polarization is 
significant (with the index of refraction differing from 
unity by more than 10%) only for B > 300-Bq, where 



B Q = m 2 e c 3 /(eh) = 4.414 x 10 13 G is the critical QED 
field strength. However, when combined with the bire- 
fringence due to the magnetized plasma, vacuum po- 
larization can greatly affect radiative transfer at much 
smaller field strengths. A "vacuum resonance" arises 
when the contributions from the plasma and vacuum 
polarization to the dielectric tensor "compensate" each 
other [H, [H, [l]]]- For a photon of energy E, the 
vacuum resonance occurs at the density p v ~ 9.64 x 
lQ~ 5 Y c r 1 B 2 2 E 2 f~ 2 g cm~ 3 , where Y e is the electron 
fraction, E\ — E/(l keV), and / = f(B) is a slowly vary- 
ing function of B and is of order unity (/ = 1 for B <C Bq 
and / -» {B/bBq) 1 / 2 for B » B Q ; see refs. [HEl). For 
p > pv (where the plasma effect dominates the dielectric 
tensor) and p < pv (where vacuum polarization domi- 
nates), the photon modes (for E <C Ese, E ^ Eb% and 
OkB 7^ 0) are almost linearly polarized; near p = pv, how- 
ever, the normal modes become circularly polarized as a 
result of the "cancellation" of the plasma and vacuum 
effects — both effects tend to make the mode linearly 
polarized, but in mutually orthogonal directions. When 
a photon propagates in an inhomogeneous medium, its 
polarization state will evolve adiabatically (i.e. following 
the K + or K- curve in Fig. ^) if the density variation 
is sufficiently gentle. Thus, a X-mode (O-mode) photon 
will be converted into the O-mode (X-mode) as it tra- 
verses the vacuum resonance, with its polarization ellipse 
rotated by 90° (Fig. QJ. This resonant mode conversion 
is analogous to the Mikheyev-Smirnov-Wolfenstein neu- 
trino oscillation that takes place in the Sun[T^. l2(ij]. For 
this conversion to be effective, the adiabatic condition 
must be satisfied [HEU (cf. pjj) 

E > E ad = 2.52 (/ tan0 fcB |l - Ul \) 2/3 ( -^J keV, 

where = (Ebi/E) 2 and H p = \dz/d\np\ is the density 
scale height (evaluated at p = pv) along the ray. For an 
ionized Hydrogen atmosphere, H p ~ 2kT/(m p g cos 0) — 
1.65 T 6 /(g u cos 9) cm, where T = 10 6 T 6 K is the tem- 
perature, g — 10 14 gi4 cm s~ 2 is the gravitational acceler- 
ation, and 9 is the angle between the ray and the surface 
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FIG. 1: The polarization ellipticity of the photon mode as 
a function of density near the vacuum resonance. The two 
curves correspond to the two different modes. In this example, 
the parameters are B = 10 13 G, E — 5keV, Y e = 1, and 
6kB = 45°. The ellipticity of a mode is specified by the ratio 
K = —iE x /E y , where E x (E y ) is the photon's electric field 
component along (perpendicular to) the k-B plane. The O- 
mode is characterized by \K\ S> I, and the X-mode \K\ <C 1. 



normal. In gene ral, the mode conversion probability is 



given by 



In gcnci 
dill 



P con = 1 - cxp [-(ir/2)(E/E ad ) 3 ] 



(2) 



The probability for a nonadiabatic "jump" is (1 — P con )- 
Because the two photon modes have vastly different 
opacities, the vacuum resonance can significantly affect 
the transfer of photons in NS atmospheres. When the 
vacuum polarization effect is neglected, the decoupling 
densities of the O-mode and X-mode photons (i.e., the 
densities of their respective photospheres, where the op- 
tical depth measured from outside is 2/3) are approxi- 
mately given by (for Hydrogen plasma and 8 k B not too 
close to 0) po ^ 0A2Tq 1/4 eI /2 G^^ 2 g cur 3 and p x ^ 
486T f r 1/4 P 1 1/2 P 14 G- 1 /2 g cm- 3 , where G = 1 - er E l kT 
[l7|. There are two different magnetic field regimes: For 
"normal" magnetic fields, 
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FIG. 2: A schematic diagram illustrating how vacuum polar- 
ization affects the polarization state of the emergent radiation 
from a magnetized NS atmosphere. This diagram applies to 
the "normal" field regime [B < 7 x I0 13 G, see eq. (3)] in 
which the vacuum resonance lies outside the photospheres of 
the two photon modes. The photosphere is defined where the 
optical depth (measured from outside) is 2/3 and is where 
the photon decouples from the matter. At low energies (such 
as E < I keV), the photon evolves nonadiabatically across 
the vacuum resonance (for 8kB not too close to 0), and thus 
the emergent radiation is dominated by the X-mode. At high 
energies (E > 4 keV), the photon evolves adiabatically, with 
its plane of polarization rotating by 90° across the vacuum 
resonance, and thus the emergent radiation is dominated by 
the O-mode. The plane of linear polarization at low energies 
is therefore perpendicular to that at high energies. 



the atmosphere structure and total spectrum can be cal- 
culated without including vacuum polarization [to be 
more accurate, we require (B/Bi) 4 <C 1 for this to be 
valid]. For concreteness, we consider emission from a hot 
spot (magnetic polar cap) on the NS; the magnetic field 
at the hot spot (with size much smaller than the stellar 
radius) is perpendicular to the stellar surface. Let the 
specific intensities of the O-mode and X-mode emerg- 
ing from their respective photospheres (which lie below 

the vacuum resonance) be Iq^ and 1% , which we cal- 
culate using our H atmosphere models developed previ- 

(0) 



B < Bi ~ 6 6 x io 13 T~ 1/8 £T 1/4 G- 1/4 G (3) ° usly ^ For a given B and TcS ' both ° and 



the vacuum resonance lies outside both photospheres 
(pv < POiPx)', for the magnetar field regime, B > Bi, 
the vacuum resonance lies between these two photo- 
spheres, i.e., po < pv < Px (the condition py < px 
is satisfied for all field strengths and relevant energies 
and temperatures). These two field regimes yield quali- 
tatively different X-ray polarization signals. 

Consider the "normal" field strengths, 10 12 G < B < 
Bi, which apply to most NSs (see Fig.[2J). In this regime, 



1^ depend on E and 9 k B at emission (the hot spot). 
As the the radiation crosses the vacuum resonance, the 
intensities of the O-mode and X-mode become Iq = 



(l-Pcon)4 0) +Pcon/r and I x = (l-P co „)/^ +P C o„/ ( 



r(0) 



(0) 



r(0) 
con J Q ■ 

In calculating P con , we use the temperature profile of the 
atmosphere model to determine the density scale height 
at the vacuum resonance. We note that in principle, cir- 
cular polarization can be produced when a photon crosses 
the vacuum resonance |2l| , but the net circular polariza- 
tion is expected to be zero when photons from a finite- 
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FIG. 3: The phase evolution of the observed spectral flux (up- 
per panel, in arbitrary units), the linear polarization Fq/Fi 
(middle panel) and the polarization degree Pl = (Fq + 

Pu) 1 ^ 2 /Fi (lower panel) produced by the magnetic polar cap 
of a rotating NS. The model parameters are: the magnetic 
field B — 10 13 G, the effective temperature of the polar cap 
Toff = 5 x 10 6 K, the angle between the spin axis and the line- 
of-sight 7 = 30° and the angle between the magnetic axis and 
spin axis f3 = 70°. The different curves are for different pho- 
ton energies as labeled. In the upper panel, the flux at 5 keV 
has been multiplied by 10 relative to the other curves. The 
definition of the Stokes parameter is such that Fq /Fi — 1 cor- 
responds to linear polarization in the plane spanned by the 
line-of-sight vector and the star's spin axis (the zero phase 
corresponds to the polar cap in the same plane). Note that 
the polarization plane of high-energy (E > 4 keV) photons 
is perpendicular to that of low-energy photons (E < 1 keV). 
Also note that for E around _E a d, the polarization degree Pl 
can be zero at certain phases such that P con = 1/2. 
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FIG. 4: The phase-averaged linear polarization as a function 
of photon energy. The model parameters are B = 10 13 G 
and T e ff = 5 x 10 6 K, and the three curves correspond to 
three different sets of angles (7, j3) as labeled. The averaged 
(Fu) = 0. The dashed curves depict the results when the 
vacuum polarization effect is turned off. 



sized polar cap are taken into account. 

To determine the observed polarization, we must con- 
sider propagation of polarized radiation in the NS mag- 
netosphere, whose dielectric property in the X-ray band 
is dominated by vacuum polarization [8|. As a photon 
(of a given mode) propagates from the polar cap through 
the magnetosphere, its polarization state evolves adia- 
batically following the varying magnetic field it experi- 
ences, up to the "polarization-limiting radius" r p i, be- 
yond which the polarization state is frozen. We set up 
a fixed coordinate system XYZ, where the Z-axis is 
along the line-of-sight and the X-axis lies in the plane 
spanned by the Z-axis and ft (the spin angular velocity 
vector). The direction of the magnetic field B following 
the photon trajectory in this coordinate system is spec- 
ified by the polar angle 9b(s) and the azimuthal angle 
4>b(s) (where s is the affine parameter along the ray). 
The X, F-components of the electric field of the photon 
mode are (Ex,Ey) = (cos 4>b, sin fis) for the O-mode 
and (— sin0s, cos 0s) for the X-mode. The adiabatic 
condition requires \nx — no\ ^ 2(hc/ E)\d<j)B /ds\ (where 
n x, no are the indices of refraction of the two modes); 
this condition breaks down at r = r p \, which is much 
greater than the stellar radius for all parameter regimes 
of interest in this paper. The observed Stokes param- 
eters (normalized to the total intensity /) are given by 
Q/I = (2P con - l)p e cos20 B (r pl ) and U/I = (2P con - 

l)p e Bin2&,(rpi), where Pe = (if - I o 0) )/(lf + I { o 0) ) 
is the "intrinsic" polarization fraction at emission, and 
4>b{j pl) is the azimuthal angle of the magnetic field at 
r = r p \ [For r p \ -C c/fi, or for spin frequency <C 70 Hz, 
^b(^pi) — 7t + 4>b(R)}- We calculate the observed spec- 
tral fluxes Fi, Fq, Fjj (associated with the intensities 
/, Q, U) using the standard procedure, including the 
effect of general relativity |2^| . 

Figure |3] shows the total flux and polarization "light 
curves" generated from the hot polar cap of a rotating 
NS. As the star rotates, the plane of linear polarization 
rotates. The light curves obviously depend on the ge- 
ometry (specified by the angles j3 and 7). Note that Fq 
for low-energy (E < 1 keV) photons is opposite to that 
for high-energy photons (E > 4 keV) , which implies that 
the planes of linear polarization at low and high ener- 
gies are orthogonal. This feature also manifests in the 
phase-averaged linear polarization (see Fig. 0J. This is 
an unique signature of photon mode conversion induced 
by vacuum polarization. 

In the magnetar field regime, B > B{, the vacuum res- 
onance lies between the photospheres of the two modes. 
At low energies (e.g., E < 1 keV), no mode conversion oc- 
curs at the vacuum resonance, and vacuum polarization 
makes the X-mode photosphere lie above (i.e. at lower 
density) its "original" location (i.e. when vacuum polar- 
ization is turned off) because the photon opacity exhibits 
a large spike at the vacuum resonance P"t1 | . At high ener- 
gies (e.g., E > 4 keV), the effective X-mode photosphere 
lies very near the vacuum resonance. For both low and 
high energies, the emergent radiation is dominated by the 
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FIG. 5: Same as Fig. 01 except for B = 5 x 10 14 G and 
T e g = 5 x 10 6 K. The different curves are for different photon 
energies (I keV dashed, 2 keV long-dashed, 5 keV solid lines). 
The thicker curves are for 7 = 30°, /3 = 70°, and the thin 
curves for 7 = 60°, /3 = 40°. In contrast to Fig. 13141 in 
the magnetar field regime, the planes of linear polarization at 
different photon energies coincide with each other. 

X-mode (for 8kB not too close to 0). Therefore the planes 
of linear polarization at different energies coincide and 
evolve "in phase" as the star rotates (see Fig. 01 . These 
polarization signals are qualitatively different from the 
"normal" field regime. Note that in the magnetar field 



regime, vacuum polarization significantly affects the total 
spectral flux from the atmosphere 0, 0, : (i) Vac- 
uum polarization makes the effective decoupling density 
of X-mode photons (which carry the bulk of the ther- 
mal energy) smaller, thereby depleting the high-energy 
tail of the spectrum and making the spectrum closer to 
black-body; (ii) vacuum polarization suppresses the pro- 
ton cyclotron line and other spectral line features [24J in 
the spectra by making the decoupling depths inside and 
outside the line similar. As we discussed in previous pa- 
pers 0, 01 > the absence of lines in the observed thermal 
spectra of several magnetar candidates |25| may be con- 
sidered as a proof of the vacuum polarization effect at 
work in these systems. Measurement of X-ray polariza- 
tion would provide an independent probe of the magnetic 
fields of these objects. 

Finally, we note that although the specific results 
shown in Figs. 13151 refer to emission from a hot polar 
cap on the NS, we expect the vacuum polarization sig- 
nature (e.g., that the planes of linear polarization at 
E < 1 keV and at E > 4 keV are orthogonal for 
B < 7 x 10 13 G) to be present in more complicated 
models (e.g. when several hot spots or the whole stel- 
lar surface contribute to the X-ray emission). This is 
because the polarization-limiting radius (due to vacuum 
polarization in the magnetosphere) lies far away from 
the star (see above), where rays originating from differ- 
ent patches of the star experience the same dipole field 
Our results therefore demonstrate the unique poten- 
tial of X-ray polarimetry 26] in probing the physics under 
extreme conditions (strong gravity and magnetic fields) 
and the nature of various forms of NSs. 
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